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Abstract. LISA is a joint space mission of the NASA and the ESA for detecting low 
frequency gravitational waves (GW) in the band IQ-^ - 0.1 Hz. The proposed mission will 
use coherent laser beams which will be exchanged between three identical spacecraft forming a 
giant (almost) equilateral triangle of side 5 x 10® kilometres. The plane of the triangle will make 
an angle of ~ 60° with the plane of the ecliptic. The spacecraft constituting LISA will be freely 
floating in the ambient gravitational field of the Sun and other celestial bodies. To achieve the 
requisite sensitivity, the spacecraft formation should remain stable, one requirement being, the 
distances between spacecraft should remain as constant as possible - that is the flexing of the 
arms should be minimal. In this paper we present a solution - the projectile solution - which 
constrains the flexing of the arms to below 5.5 metres/sec in a three year mission period. This 
solution is obtained in the field of the Sun and Earth only, which principally affect the motion 
of the spacecraft, especially the flexing of LISA'S arms. 
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1. Introduction 

LISA - Laser Interferometric Space Antenna - is a proposed mission of the ESA and NASA 
which will use coherent laser beams exchanged between three identical spacecraft forming a 
giant (almost) equilateral triangle of side 5 x 10^ kilometres for observing low frequency cosmic 
GW [l]. This will complement the ground-based detectors which are geared to operate at higher 
frequencies ranging from few tens of Hz to kHz. For the successful operation of LISA it is crucial 
that the formation of spacecraft be stable - that is, the spacecraft should maintain as much as 
possible, constant distances between them. However, the spacecraft are freely floating in the 
ambient gravitational field of the Sun, planets and other celestial bodies (moon for instance) and 
it is a astrometry problem to seek spacecraft orbits which maintain the equilateral triangular 
formation as nearly as is possible - that is, optimal orbits for the spacecraft should be found. 
There are several criteria which the spacecraft formation should satisfy for LISA's successful 
operation - constraints on, variation in armlengths, the angles between arms, etc. Here we focus 
on the variation in armlengths, the so-called 'flexing' of the arms for the reasons detailed below. 
Optimisation of LISA orbits will be also useful in simplifying the hardware that will be required 
in the design of LISA. 

Minimising the flexing of the arms is important for suppressing the laser frequency noise. 
In ground-based detectors, the near exact symmetry between the arms suppresses this noise as 
it is common to both arms. But in LISA such high symmetry is not possible, and moreover, 
the armlengths change with time. Suppression of this noise is crucial since the raw laser noise 
is orders of magnitude larger than other noises in the interferometer. In LISA, six data streams 
arise from the exchange of laser beams between the three spacecraft. The cancellation of the noise 
is achieved by the technique called time-delay interferometry (TDI) where the six data streams 
are combined with appropriate time-delays [2]. This is possible because of the redundancy 
present in the data. TDI was put on a sound mathematical footing by establishing that the data 
combinations had an algebraic structure. The time delayed data is represented by polynomials 
of time-delay operators acting on the data, each time-delay operator playing the part of an 
'indeterminate' of a polynomial ring. The data combinations are then represented by polynomial 
vectors which form a free module over the polynomial ring of time-delay operators. Out of 
these, the data combinations cancelling laser frequency noise form a submodule of this free 
module, wellknown in mathematics, as the first module of syzygies j^. The generators of the 
submodule were found assuming constant armlengths, where one then deals with the simpler 
case of a commutative ring of time-delay operators [Sj |4]. But for realistic spacecraft orbits, 
the armlengths vary with time, and then the TDI methods involve non-commutative operators 
leading to the imperfect cancellation of laser frequency noise or the presence of residual noise. 
The residual noise in turn depends on the rate of change of armlengths - the flexing of arms; 
thus searching for orbits which reduce the flexing also reduces the residual laser frequency noise. 

In this paper, following [5], we include the gravitational fleld of the Earth in addition to 
that of the Sun's in the optimisation problem. The orbits in the Sun's fleld are taken upto 
second order in a (or eccentricity), where a = 1/2R, where Z '--^ 5 x 10^ km is the nominal 
distance between the spacecraft and R is one astronomical unit ~ 1.5 x 10^ km. We flnd the 
perturbative approach for the Sun's fleld convenient because we also introduce the Earth's effect 
perturbatively. The second order terms in a involve the Sun's fleld upto the octupole order 
and as shown in [6] almost exactly replicate the Keplerian orbits of the spacecraft and therefore 
also the flexing. We then linearly superpose the perturbative effect of the Earth's gravitational 
fleld over the Sun's fleld. We choose the Earth over .lupiter because the Earth perturbs the 
Keplerian orbit in resonance, resulting in unbounded growing of the perturbations and also as 
shown in [5j Jupiter's tidal fleld which affects the flexing is less than 10% of the Earth's and hence 
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not a dominant one. Although we recognise that the problem is inherently non-linear - it is a 
three body problem - the linear perturbative approach we believe will be useful for short mission 
periods and also provide directions towards solving the fully general optimisation problem. The 
analytic approach which we follow here helps to gain insight into the problem. 

We first obtain the general solution containing 18 arbitrary constants, corresponding to 3 
positions and 3 velocities for each of the three spacecraft. Optimising the 18 parameter solution is 
a daunting problem - we do not attempt to do so here. However, from physical considerations, we 
present a solution, which we call the 'projectile' solution which considerably reduces the flexing 
of the arms - the rate of change of all armlengths is less than 5.5 metres/sec in a three year 
mission period. We believe that these insights will lead us to the full solution of optimisation on 
18 parameters. (The nomenclature 'projectile' solution will be justified later in the text). 

2. The general perturbative solution of the spacecraft orbits 

In the subsection below, we briefly summarise the results of the previous papers [5, 6j and write 
down the perturbed CW equations in the two small parameters which describe the effects of the 
Sun and Earth. In the next subsection we write down the general solution with 18 arbitrary 
constants; these are the 18 parameters to be varied in order to optimise LISA's orbits with 
respect to given criteria. Here, our sole criterion is minimising the flexing of the arms. 

Alternatively, we could have chosen to work with the exact Keplerian orbits of the spacecraft 
for the Sun's field and added to these the perturbations due to the Earth as we have done in [5]. 
However, here we use the approximate solution to the second order in the parameter a because 
with this solution we gain important physical insights into the problem - for example, we can 
slightly adjust the tilt of the plane of LISA, by choosing certain constants judiciously, which then 
helps in reducing the flexing of the arms. Moreover, we have shown in [6], we lose very little in 
accuracy, because the approximate solution is extremely close to the exact. Also it is gratifying to 
check that, although we end up with 12 arbitrary constants for each spacecraft in this approach, 
they combine two by two, to yield only six independent arbitrary constants corresponding to the 
six initial conditions on the three position coordinates and three components of the velocity. 

2.1. The perturbed Clohessy- Wiltshire equations 

If we consider only the Sun's field, there exist orbits in which the plane of the LISA triangle 
makes an angle of about 60° with the ecliptic and the cluster rolls once per year and for which the 
armlengths remain constant upto a percent. For these orbits, to the first order in the eccentricity, 
the distances between spacecraft remain constant; only at the second order in eccentricity the 
variations in armlengths appear. It was shown in this case that the flexing could be reduced to a 
minimum ~ 48, 000 km j6] by judiciously choosing the orbital parameters of the spacecraft. For 
establishing this result, it was found convenient to use the Clohessy- Wiltshire equations [?]• 

Clohessy and Wiltshire make a transformation to a frame - the CW frame {x, y, z} which 
has its origin on the reference orbit and also rotates with angular velocity The a; direction 
is normal and coplanar with the reference orbit, the y direction is tangential and comoving, 
and the z direction is chosen orthogonal to the orbital plane. They write down the linearised 
dynamical equations for test-particles in the neighbourhood of a reference particle (such as the 
Earth). The length scale here is the Earth-Sun distance of 1 A. U. and the motion of a test 
particle is described by these equations if its distance from the origin is small compared with this 
length scale. Since the frame is noninertial, Coriolis and centrifugal forces appear in addition to 
the tidal forces. 
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We take the reference particle to be orbiting in a circle of radius R with constant angular 
velocity fl. Then the transformation to the CW frame {x,y,z} from the barycentric frame 
{X, Y, Z} is given by, 

X ^ {X - Rcos nt) cos nt + {Y ~ R sin fit) sin fit , 
y = - (X - Rcos nt) sin fit + {Y ~ R sin ilt) cos fit , 

z ^ Z. (1) 
The unperturbed CW equations for a test particle with coordinates (x, y, z) are given by, 
i-2ny- sn^x = , 
y + 2fii = , 

z + n'^z =0. (2) 

These equations include terms upto the quadrupole, when the Sun's field is Taylor exanded about 
the origin of the CW frame. The solutions to these equations we call the zero'th order. Among 
these we choose the solutions which form an equilateral triangular configuration of side I. For 
the fcth spacecraft. A: = 1, 2, 3 we have the following coordinates: 

Xk^ - ipo cos{nt - 2n{k - l)/3 - ^o) , 
yk = Po sin{flt - 2TT{k - l)/3 - ^o) , 

Zk = ^po cos(m - 27r(fc - l)/3 - ^o) , (3) 

where po = l/^/i is the constant distance each spacecraft maintains from the origin of the CW 
frame and is an arbitrary constant phase. In this solution, any pair of spacecraft maintain 
the constant distance I between each other. 

In we have shown that if we include the octupolar terms and solve perturbatively using 
the zeroth order solution as given by Eq.|(3]), we obtain the fiexing of the arms due to the Sun's 
field only. We now include the Earth's field as well. LISA follows the Earth 20° behind. We 
consider the model where the centre of the Earth leads the origin of the CW frame by 20° - 
thus in our model, the 'Earth' or the centre of force representing the Earth, follows the circular 
reference orbit of radius 1 A. U. Also the Earth is at a fixed position vector = (x©, y®, z^) in 
the CW frame. We find that a;® = -R{1 - cos 20°) - -9 x 10^ km, = i?sin20° - 5.13 x 10^ 
km and = 0. In order to write the CW equations in a convenient form we first define the small 
parameter e in terms of the quantity — GM^/d^, where c?® = jr®! ~ 5.2 x 10^ km is the 
distance of the Earth from the origin of the CW frame; we define e — >~^%l^ ~ 7.16 x 10~^ which 
is essentially the ratio of the tidal force exerted by the Earth to that of the Sun. We approximate 
|r — rg| by d® in the force field of the Earth. We then linearly add the two perturbative terms, 
namely, the terms describing the octupolar field of Sun and the Earth's field and obtain the 
perturbed CW equations: 

i-2VLy~- i^fx + 5^(2x2 -y^ - z^) + e^^ix - x®) = , 
y + 2nx J — xy + ef]^(y - ?/®) =0, 

z + n^z J — xz + efl'^z = 0. (4) 

We now have the perturbed equations in two small parameters a and e. We seek perturbative 
solutions to Eq. ^ to the first order in a and e. We note that the forcing terms by the Earth 
in these equations appear at the same frequency fl and hence they imply resonance. This means 
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that the Earth's effect on LISA is cumulative and therefore important. Also, we have here ignored 
higher order terms in both a and e as well as the cross terms in these parameters, in order that 
the problem becomes linear and therefore tractable. Thus the solutions we will obtain are valid 
in the short term or for short periods of the LISA mission; for longer periods the problem is 
inherently nonlinear and difficult to deal with analytically - it is infact a three body problem 
that we are approximating. 

The quantity lu^^ ~ 18.8 years defines a timescale. If we assume a stationary Earth and 
Earth's gravitational field only, the free fall time of a particle initially at rest at a distance 
from Earth is about 21 years which is comparable to this timescale. The three year mission 
period we have assumed here is smaller than the above timescales and therefore sufficiently short 
for our analysis to be useful. Moreover, as it will turn out for the solution we present, the 
LISA spacecraft fall towards the Earth about half a million km from their initial positions, thus 
remaining well within the CW frame. Thus we expect the linear perturbative analysis that we 
have carried out here to hold good. 



2.2. The general solution with 18 parameters (arbitrary constants) 

Since these solutions have been derived in previous papers [5][^, we merely state the results 
here. We adopt the following notation: we denote the solution for spacecraft k with k = 1,2,3 
by the bracketed suffix k and the zeroth order by the suffix 0, the a perturbation by the suffix 
1 and Earth's perturbation - the e perturbation - by the suffix 2; thus we write: 

X(k)=X(^k)^ + ax(k)i+^X(k)2, (5) 

and similarly for the y and z coordinates. These are the coordinates of the spacecraft in the CW 
frame. Further, to reduce the clutter, we choose units of time and length such that SI = 1 and 
1 = 1. In these units, to the zeroth order, the spacecraft form an equaliteral triangle of side unity 
with the distance of each spacecraft from the origin equal to 1/ \/3; also one year period in these 
units equals t = 27r. In these units we may rewrite Eq.(l3| as follows: 



1 

2 

1 

1 



2/(fc)o = sm (pk ■ 



^{k)0 



2 cos (j3k , (6) 



where, 4>k =t — 27r(fc — l)/3 — to and t^ is a arbitrary constant phase. The perturbative solutions 
are the following: 

5 1 

x^k)i = 2ylfe + Bk cos 0fc + Cfe sin 0fc + - - — cos 20^ , 

o 24 

V(k)i = - {'iAk + 5/4) t + 2{Ck cos (l>k - Bk sin (l)k) + Dk + ^ sin 2(j)k , 

\/3 1 

Z{k)i = Ek cos (pk + Fk sin 4>k + ^ cos 20^ ; (7) 



and, 



C(fe)2 = 2A'f. + X® + 2%© + B'^ cos (j)k + C'f, sin (f)k + sin (j)k 



5t 
475' 



2/(fe)2 = -3At- 2a;©i - -y®t + cos 0fe - — sm (j)t 
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Z(k)2 = E'k COS 0fc + Fl, sin (f)k + -t sin 0^ . (8) 

The quantities Ak, Bk,Ck, Dk, Ek, and A',., B'i^,C'f., D'/,, E'l^, Fj. are arbitrary constants. For 
each spacecraft, there seem to be 12 arbitrary constants. However, if we now add up all the 
solutions given in Eqs.(l6]), ((?]) and ^ to obtain the full solutions as in Eq.|[5]), the arbitrary 
constants combine as aAk + eA'j^,aBk + eB'f,,... etc. to give just six independent arbitrary 
constants for each spacecraft as demanded by the three second order simultaneous differential 
equations. We find however, that it is better to leave the arbitrary constants as they are, 
because from our previous experience, we know what values the arbitrary constants should take 
in order that the spacecraft form stable or nearly stable configurations in which the variation in 
armlengths is acceptably small. 

3. Stability and reduced flexing of the arms of LISA: the projectile solution 

We seek solutions that are (i) stable, and (ii) reduce the fiexing of the arms. Following [6] we 
satisfy the first criterion by choosing the following values of the constants: 

Ak = - — , Bk ^—,Ck = Dk = 0, Ek = ^, Fk^O, k = 1,2,3. (9) 
12 Id Id 

This choice of constants ensures that, (i) to the first order in eccentricity (or a), the spacecraft 

maintain constant distances from the origin, forming an equilateral triangle which makes an angle 

of 60° with the ecliptic; (ii) to the second order in a, the spacecraft do not drift away - the choice 

of Ak ensures that the secular term proportional to t in the y(k)i is set to zero; now the angle of 

the plane of the triangle is not exactly 60°, but very close to it - 5a/16 < 0.01 radians from 60°. 

If we 'switch off' the Earth's field, this choice of tilt angle (of constants) ensures that the fiexing is 

kept at a minimum to < 48, 000 km which is less than 1% variation in the armlength. Thus, even 

with the Earth's field we seek solutions that are close to the previously found solutions, which 

were shown to have optimal properties in the field of the Sun only. With this 'safe' strategy, we 

expect not to stray away from optimality. The solutions that we will find do not lay any claim 

to exact optimality, but they do exhibit adequate reduced fiexing of the arms to < 5.5 m/sec in 

a 3 year mission, which would inturn be useful for reducing the residual laser frequency noise in 

TDI. In [5] it has been shown that in the type of solutions we are considering, it is essentially 

the L terms which contribute to the residual noise, where L generically is the length of any one 

arm; higher order terms can be neglected. Moreover, the amplitude of the residual noise is cx i 

and thus its power spectral density (PSD) is cx L^. Thus a reduction in L from 10 m/sec say, 

which was the estimate in earlier literature, to 5.5 m/sec, reduces the PSD of the residual noise 

to almost 30% of its earlier estimate. 

We now turn to the primed set of constants which occur in the Earth's perturbative part 

of the solution. In this part we are guided by a physical criterion. We would Hke the Earth's 

perturbative effect to be small during LISA's mission. One way to achieve this is by setting 

X{k)2 = y(k)2 = ^(fe)2 = as also the velocities i;(/c)2 = y{k)2 = ^(k)2 = at an appropriate epoch 

t, where the overdot represents the time derivative of a quantity. This appropriate epoch we 

choose at the middle of the mission; so if the mission period is T and we arrange so that the 

appropriate epoch occurs at t = 0, then the mission duration is —T/2 < t < T/2. With these 

initial conditions the primed set of arbitrary constants are determined. Defining for the spacecraft 

A: = 1, 2, 3, the constant phases tfe to + 27r(fc — l)/3, and imposing the above mentioned initial 

conditions at t — 0, the arbitrary constants for the spacecraft take the values: 

1 
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2 

73 



, cos tk - 2y(5 sin tk 



V3 



sin tk 



C'k = X© sintk - 2y( 



V3 . 

, cos tk — sm tk cos tk ; 



D'k = 42/g 



E'k 



F'k 



-Isin^,, 

tk , 



— sm 

4 

1 

- sm tk cos tk ■ 

4 



(10) 



These constants determine the orbits of the spacecraft and therefore the distances between them 
as a function of the epoch t. The armlengths kj, i, j — 1,2,3 are then computed by the formulae: 



= - ^{J)f + iVi'i) - y(j)f + - Z(j)f 



(11) 



In the figures |T]) and ^ below, we plot the armlengths Uj and the rate of change of armlengths 
lij for the phase to = for the mission period of 3 years which in our units is T = 67r. The 
armlength variation increases from the previous optimum obtained in the field of the Sun only, 
from 48,000 km to roughly 60,000 km in the combined field of the Sun and Earth. Also the rate 
of change of armlengths increases from the maximum of 4 metres/sec for Sun's field only to a 
maximum of 5.5 metres/sec in the combined field. 




Figure 1. The figure shows the variation in the three armlengths of LISA for a mission 
period of three years (— Stt < fit < Stt) for the phase to = in millions of km. The maximum 
variation in armlengths is about 60,000 km. 

We also plot in figure ([3]) the rate of change of armlength Z12 for the different phases, 
to = 0,40°, 80°. We observe that the profiles of the curves essentially repeat; the curves are 
basically time translated. The figure shows that the maximum fiexing is essentially insensitive 
to the phase. This is the consequence of the high symmetry of the LISA configuration. We 
term this as the 'projectile' solution because just as when a stone (projectile) is vertically thrown 
from the ground in a gravitational field that is assumed to be constant, the stone reaches zero 
velocity midway and returns to the ground. The solution presented here describes an analogous 
situation; LISA is 'thrown' away from the Earth initially, thus it moves away relative to the 
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Figure 2. The figure shows the time derivative of the armlengths for a mission period of 
three years (— Stt < Qt < Stt) for the phase to = 0. The flexing of the arms is less than 5.5 
metres/sec. 



Phase 0° 
Phase 40° 
Phase 80° 




Figure 3. The figure shows the flexing of the arm h2 for the different phases to = 0,40°, 80° 
and — Stt < fit < 3tt. The maximum flexing is seen to be insensitive to the phase. 



Earth in the first half of the mission, reaches maximum distance away from the Earth at mid- 
mission, and then in the next half of the mission period falls towards the Earth. One can easily 
compute the distance travelled by LISA relative to the Earth by examining the expressions 
for {x2,y2,Z2) in Eq.®. We have dropped the subscript k in order to avoid clutter - we can 
do this because of symmetry of the LISA configuration; the results are essentially the same 
for all spacecraft. Let us therefore consider spacecraft 1. If we consider a 3 year mission 
period, then T — Gtt. At initial time and final times, t — ±T/2 — ±37r, the coordinate which 
dominates is the y coordinate (not surprisingly, as this is roughly the direction of the force of 
the Earth); and the term that dominates in the y coordinate is the quadratic term in t. Thus 
y2(±r/2) - (3/2?/e)(rV4) = 277rV22/e ^ 1367 in the units chosen, for T = 67r. Converting 
to km by multiplying by the factor e x 5 x 10^ km yields ey2 ~ 4.9 x 10^ km. At t = 0, the 
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initial condition implies, 2/2 = and also 1/2 = 0. Thus in this solution, in the first year and 
half, LISA travels about 500,000 km away from the Earth and then falls back towards the Earth 
about the same distance in the second half of the mission. It does not fall back exactly to the 
same point though, even relative to the Earth, because X2 ~ ^y^t for large times and thus 
ea;2(±37r) ~ ztGney^ ~ ±7 x lO** km. The z coordinate changes very little and plays a minor 
role in the solution. 

As compared to the solution described in [5J, in which the initial conditions a;(fc)2 = 2/(/c)2 — 
Z{k)2 — i(k)2 — y(k)2 = Z(k)2 = wcrc applied at the start of the mission t — 0, and where it was 
found that the maximum flexing after three years was about 8 metres/sec, here in the projectile 
solution, the flexing is reduced to 5.5 metres/sec, in which the same initial conditions are applied 
mid-mission. This gives an improvement of about 30% in the maximum flexing, and more than 
a factor of two in the PSD of the residual laser frequency noise. If this level of residual noise 
in the TDI observables can be tolerated, then these spacecraft orbits can be considered to be 
adequate. 

4. Concluding Remarks 

In this paper, we have presented a solution for the LISA spacecraft orbits which gives reduced 
flexing of LISA's arms to less than 5.5 metres/sec in a three year mission. The solution has 
been obtained in the combined field of Sun and Earth. The solution although approximate is 
analytical and hence provides valuable insights into the problem of optimisation. Clearly, this is 
not the most optimised solution that is possible. The truly optimised solution for short mission 
periods may be computed by varying the 18 arbitrary constants in the general solution which 
has been given in this paper. Although we have argued here, that the tidal effects due to Jupiter 
are small, for a complete solution, it would be desirable to include the field of Jupiter in future 
endeavours. 

Optimisation of orbits is an important problem for LISA, because the judicious choice of 
orbits can lead to several advantages. As we have argued here, reducing the flexing of the arms 
from say 10 metres/sec to 5.5 metres/sec, tends to reduce the PSD of the residual noise in the 
TDI observables to about 30% of its original estimated value. This reduction could further 
help in the simplification of the TDI in which the first generation modified TDI could suffice, 
thus in turn reducing the degree of the polynomials in the time delay operators. Lower degree 
polynomials are preferred because they decrease the interpolations required to be carried out on 
the data and in turn the overall noise. Further, the optimisation of orbits can also help in the 
simplification of hardware in the design of LISA. 
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